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A PERTURBATION  ANALYSIS  OP  THE  ATTENUATION  AND 
DISPERSION  OF  SURFACE  WAVES 

H.F.  Tiersten  and  B.K.  Sinha 
Department  of  Mechanical  Engineering, 
Aeronautical  Engineering  & Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 

ABSTRACT 


A perturbation  formulation  of  the  equations  of  linear  piezoelectricity  is 
obtained  using  a Green's  function  approach.  Although  the  resulting  equation 
for  the  first  perturbation  of  the  eigenvalue  strictly  holds  for  real  perturba- 
tions of  real  eigenvalues  only,  it  is  formally  extended  to  the  case  of  purely 
imaginary  perturbations  of  real  eigenvalues.  The  extended  equation  is  applied 
in  the  calculation  of  the  attenuation  of  surface  waves  due  to  the  finite  elec- 
trical conductivity  of  thin  metal  films  plated  on  the  surface  and  air  loading. 
The  influence  of  the  viscosity  of  the  air  is  included  in  the  air  loading 
analysis,  and  the  calculated  attenuation  increases  accordingly.  Since  the 
metal  films  are  thin  compared  with  a wavelength,  an  approximate  thin  plate 
conductivity  equation  is  employed  in  the  determination  of  the  attenuation  due 
to  the  electrical  conductivity  of  the  films.  The  resulting  attenuation  is 
obtained  over  a very  large  range  of  values  of  sheet  conductivity.  This  is 
accomplished  by  using  the  equation  for  the  first  perturbation  of  the  eigenvalue 
iteratively  to  determine  the  solution  and  attendant  attenuation  to  any  desired 
degree  of  accuracy.  The  phase  velocity  dispersion  curve  due  to  the  mechanical 
effect  of  a thin  film  plated  on  a substrate  is  determined  for  relatively  large 
wavelengths  by  employing  the  perturbation  equation  iteratively,  and  excellent 
agreement  is  obtained  with  the  results  of  other  more  direct  approaches.  The 
calculations  have  been  performed  for  an  aluminum  film  on  either  ST-cut  quartz 
or  Y-Z  lithium  niobate. 


1.  Introduction 

The  attenuation  of  acoustic  surface  waves  propagating  along  electroded 
piezoelectric  substrates  arises  from  a number  of  different  causes.  In  a manner 
of  speaking,  these  can  be  loosely  categorized  as  substrate  and  electrode  mater- 
ial attenuation  and  air  loading.  The  material  attenuation  is  a result  of 
material  viscosity  and  electrical  and  thermal  conductivity.  In  the  substrate 
the  viscosity  is  most  important,  the  thermal  conductivity  of  secondary  import- 
ance and  the  electrical  conductivity  is  negligible.  In  the  thin  metal  film  the 
electrical  conductivity  is  of  primary  importance  and  the  viscosity  and  thermal 
conductivity  are  of  secondary  importance.  In  any  event  in  this  work  the  atten- 
uation due  to  the  finite  electrical  conductivity  of  the  thin  electrode  plating 

and  air  loading  are  determined.  Measurements  of  attenuation  due  to  a number 

1-5 

of  the  aforementioned  causes  have  been  made 

The  attenuation  of  surface  waves  due  to  air  loading  was  first  treated 
analytically  by  Campbell  and  Jones^,  who  ignored  the  influence  of  the  viscosity 
of  the  air  and,  consequently,  obtained  values  of  attenuation  below  those 
measured  by  Slobodnik^.  Subsequently,  Auld7  obtained  results  identical  with 
those  of  Ref. 6 using  his  first  order  piezoelectric  perturbation  theory.  The 
attenuation  of  acoustic  surface  waves  due  to  an  adjacent  semiconducting  film 

g 

has  been  calculated  by  Ingebrigtsen  by  means  of  a perturbation  procedure. 

In  this  paper  a piezoelectric  perturbation  theory  is  obtained  from  a 

9 10 

Green's  function  formulation  of  the  equations  of  linear  piezoelectricity 

The  resulting  equation  for  the  first  perturbation  of  the  eigenvalue  is  applied 
in  the  determination  of  the  attenuation  of  surface  waves  due  to  the  finite 
electrical  conductivity  of  thin  metal  films  plated  on  the  surface  and  air  load- 
ing. The  directly  determined  time  attenuation  is  converted  to  spatial  attenua- 
tion. The  resulting  perturbation  expression  for  the  spatial  attenuation  is 


2. 


similar  to  that  of  Auld7.  In  the  case  of  attenuation  due  to  air  loading  the 
stress  in  the  air  due  to  the  unperturbed  surface  wave  displacement  field  is 
inserted  in  the  equation  for  the  first  perturbation  of  the  eigenvalue.  The 
influence  of  viscosity  is  included  in  addition  to  fluid  pressure  in  the  analysis 
of  the  air.  As  a consequence  the  calculated  attenuation  due  to  air  loading  is 
larger  than  that  of  either  Ref. 6 or  7 and  is  in  excellent  agreement  with  the 
measurements  of  Szabo  and  Slobodnik^  for  ST-cut  quartz  and  Slobodnik'’  for 
lithium  niobate. 

Since  the  aforementioned  metal  films  are  thin  compared  with  a wavelength, 

12 

an  approximate  thin  plate  conductivity  equation  is  obtained,  which  enables 
the  entire  electrical  conductivity  effect  of  the  plating  to  be  treated  as  a 
boundary  condition  at  the  surface  of  the  substrate.  This  approximate  thin  plate 
conductivity  equation,  but  with  a sign  difference,  has  been  employed  by  Adler^ 
in  a treatment  of  semiconducting  thin  films.  The  attenuation  resulting  from 
the  conductivity  of  the  thin  film  is  determined  by  employing  the  aforementioned 
thin  plate  conductivity  equation  in  the  equation  for  the  first  perturbation  of 
the  eigenvalue  in  the  usual  manner.  However,  in  this  attenuation  analysis  the 
perturbation  equation  is  used  successively  to  increment  the  solution  at  one 
value  of  sheet  conductivity  from  the  previous  solution  at  a nearby  value  of 


sheet  conductivity,  because  we  are  interested  in  the  attenuation  over  a large 
range  of  values  of  sheet  conductivity.  In  addition,  the  perturbation  equation 
is  used  iteratively  at  a fixed  value  of  sheet  conductivity  to  determine  the 
solution  and  attendant  attenuation  to  any  desired  degree  of  accuracy.  With 
this  latter  procedure  small  changes  in  the  phase  velocity  are  obtained  as  well 
as  the  attenuation.  However,  we  find  that  in  this  manner  we  are  unable  to 


p 


r 


3. 


arises  because  in  that  range  the  attenuation  is  too  large  for  the  perturbation 
procedure  we  employ. 

In  addition  to  the  foregoing  the  perturbation  equation  is  applied  in  the 

determination  of  the  initial  slope  of  the  phase  velocity  dispersion  curve  due 

to  the  mechanical  effect  of  the  thin  film  and  excellent  agreement  is  obtained 
14 

with  other  more  direct  approaches.  Again,  the  iterative  procedure  is  employed 
to  obtain  considerably  more  than  the  initial  slope  of  the  dispersion  curve  from 
the  perturbation  equation,  and  the  results  are  shown  to  be  extremely  accurate. 
Finally,  it  should  be  noted  that  the  perturbation  theory  is  applicable  in  the 
determination  of  the  attenuation  due  to  all  the  aforementioned  other  causes. 

In  particular,  the  attenuation  resulting  from  the  material  viscosity  of  the 
substrate  can  readily  be  evaluated  in  the  same  manner  using  the  results  of 
Lamb  and  Richter^  for  quartz. 

2 . piezoelectric  Perturbation  Theory 

In  this  section  we  obtain  a piezoelectric  perturbation  theory  from  a 
Green's  function  formulation''’*”'  of  linear  piezoelectricity.  To  this  end  we 
write  the  equations  of  linear  piezoelectricity  in  the  form1*’ 


T.  . . + T.  . . = pu  . , 

13,1  i3,i  3 

D 


. + D.  . = 0 , 
1,1  1,1 


(2.1) 

(2.2) 


i 


Tij  Cijk£Uk,-C  + ekij^,  k’ 


Di  " eik£Uk,  l ~ ‘ik'P.k 


(2.3) 


where  (2.1)  and  (2.2)  are  the  stress  equations  of  motion  and  charge  equation 
of  electrostatics,  respectively;  and  (2.3)  are  the  linear  piezoelectric  con- 
stitutive equations  and  T„,  Uj  and  th  denote  the  components  of  stress,  mech- 
anical displacement  and  electric  displacement,  respectively;  p and  cp  denote 


L 


■: 

L. 


4. 


the  mass  density  and  electric  potential,  respectively;  c anc*  ejj 

denote  the  elastic,  piezoelectric  and  dielectric  constants,  respectively. 

G G 

The  quantities  and  denote  volumetric  mechanical  and  electrical  perturba- 
tion terms,  respectively.  In  (2.1)-  (2.3)  we  have  employed  the  convention  that 
a comma  followed  by  an  index  denotes  partial  differentiation  with  respect  to  a 
space  coordinate,  the  dot  notation  for  differentiation  with  respect  to  time  and 

the  summation  convention  for  repeated  tensor  indices. 

17 

The  equations  for  the  Green's  tensor  corresponding  to  (2.1)  - (2.3)  may 


be  written  in  the  form 


a 2a.  .a 

T.  . . + puj  g.  = - 6 (p-  Q)6  . , 
13,1  3 3 ’ 

<i  =-6<p-K> 

JX  E c01 

T.  . = c.  .G,  . + e,  . .f  , , 

ij  rjkX  k,  l kig  , k ’ 

a s ,a 

i eik^Gk,i  " eikf,k  ’ 


where  a,  3 = 1-4,  P and  Q denote  the  fixed  field  point  and  variable  source 

Ot 

point,  respectively,  6 is  the  Dirac  delta  function,  6^  is  the  Kronecker  delta, 

Q(  Oi 

Gj  and  f are  the  mechanical  displacement  Green's  tensor  and  electric  potential 
Green's  function  plus  cross  terms,  respectively,  and  the  significance  of  the 
remaining  quantities  in  (2.4)-  (2.6)  is  obvious  from  (2.1)-  (2.3).  In  (2.4)  - 
(2.6)  we  have  assumed  that  all  variables  have  a time  dependence  of  e1CUt.  We 
now  make  the  same  assumption  in  (2.1)  - (2.3)  and  in  the  usual  manner,  from 
(2.1)  and  (2.4),  we  form 

f T (T  . +T®  +pu)2u  )G™~  (T™  +pu)2G%  6(P-Q)6™)U  1 dV(Q)  =0,  (2.7) 

JJ  J JJJ 


kJH 


where  m = 1-3  and  all  variables  in  (2.7)  have  spatial  dependence  only. 
Employing  the  divergence  theorem  and  using  Eqs.  (2.2),  (2.3),  (2.5)  and  (2.6), 

we  obtain 


5. 


u (P)  = [ n.l  T.  .G1"  - Tm.u . + D.f™  - dS(Q)  + 

m .1  1L  13  3 13  3 1 1 J 

S 

f f T°  . ,Gm  + D6  .r1"]  dV(Q)  . 

J L 13,1  3 1,1  J 

Similarly,  from  (2.2)  and  (2.5)  we  form 

f f (D . . + D?  ,)f4-  (D4  . + 6 (P  - Q)  )cp  1 dV(Q)  = 0 , 

3 L i,  1 1, 1 1,  a — J 


(2.8) 


(2.9) 


and  utilizing  the  divergence  theorem  and  (2.1),  (2.3),  (2.4)  and  (2.6),  we  obtain 

r r 4 „4  441 

(P)  = J n^Djf  - -fyp  + T^G.  - T u J dS(Q)  + 


cpl 


JKi 

V 


.4  e 4' 
f + T. . .G. 
13,1  3 


dV  (Q ) . 


(2.10) 


Equations  (2.8)  and  (2.10)  constitute  the  Green's  function  (or  tensor)  form  of 
the  equations  of  linear  piezoelectricity.  It  turns  out  that  in  the  sequel, 
although  we  have  great  use  for  (2.8)  we  have  no  use  for  (2.10)  because  of  the 
particular  type  of  perturbation  problem  of  interest  here. 

Vibrational  eigensolutions  of  the  equations  of  linear  piezoelectricity 
satisfy  an  equation  which  may  be  written  in  the  form 


(*  n.  T ,uV  - TV  .u^  + D^CpV  - DVcp^l  dS  = (u)2-u)2)  [*  pul‘.luv  dV  , (2.11) 

J iL  13  3 13  3 iT  V J v M-  J 3 3 


where  the  scripts  p and  -g  refer  to  the  Pth  and  \>th  eigensolutions,  respectively. 
It  should  be  noted  that  (2.11)  holds  even  if  intermediate  surfaces  of  discon- 
tinuity exist.  Clearly,  from  (2.11),  for  homogeneous  boundary  conditions,  we 
have 

V dV  = N~  .6  , 

3 (II)  PV  ’ 


J pUjU 


(2.12) 


6. 


which  is  the  orthogonality  condition  for  piezoelectric  vibrations.  We  now  assume 

p n, 

that  a complete  set  of  eigensolutions  u_.  and  cp  exists  and  define  orthonormal 
eigensolutions  to  our  eigenvibration  or  eigenwave  problem  by 


j (lO 


(2.13) 


We  now  expand  the  mechanical  displacement  Green's  tensor  and  electric 
potential  Green's  vector  fm  in  the  forms 


Gm  = y M'yt . fm  = y.  Mm  & 

3 L I*  ] U M- 

p p, 


(2.14) 


where  g^  and  7"  constitute  orthonormal  solution  functions  satisfying  the  approp- 
riate homogeneous  form  of  (2.1)  - (2.3)  subject  to  the  appropriate  boundary 

conditions.  Substituting  from  (2.14)  into  (2.4),  employing  (2.6)  and  the  homo- 

v 

geneous  form  of  (2.4)  for  every  |i,  contracting  with  g^,  integrating  over  V and 
utilizing  (2.12),  we  obtain 


= - w2)  . 

p m p. 

Substituting  from  (2.14)  and  (2.15)  into  (2.8),  we  obtain 


(2.15) 


» -V 


9m(P) 


m L .2  2 ! J i i j j 

jj,  (10^-it  ) s 


F [ n (T.  / (2)  - u if  (Q)  + D.?  (Q)  - (8)  1 dS  (Q) 

u>JiiJ3  3 13  1 1 


+ I [T®.  • (Q ) + D®  r (Q)ldV(Q)  . 
J 13 , 1 3 1,  1 


(2.16) 


We  may  now  obtain  our  perturbation  procedure  from  (2.16)  in  the  usual  way  , 
i.e,,  by  letting  um  (p)  be  very  near  one  of  the  g^,  say  Then  we  may  write 


Um  = gm(P)  + 7 gi<P)V%2‘U,2>  ’ 

6.71 


(2.17) 


H = [ n . Ft  . . c/t  - u . . + D.f^i-  cpj^tldS  + j Ft?  . .g'l'  + .?L  ] dv  , (2. 18) 

p.  J iL  13^3  3 13  1 1J  J L 13,1  3 1,1  J 


1 


7. 


and  we  have 


Hl/  (lUl  ~ u,2)  = 1 > 


(2.19) 


which  is  the  equation  for  the  first  perturbation  in  eigenf requency.  If 


A = CD, 


(2.20) 


from  (2.19),  we  have 


A = H1/2o)1 


(2.21) 


for  the  first  perturbation  in  eigenfrequency.  If  we  let  the  phase  velocity  v 
be  written 


v = v^  - e 


(2.22) 


then  for  constant  wavenumber  § = § from  (2.20)  and  (2.22),  we  have 

e = A/§1  , 

for  the  perturbation  in  phase  velocity.  Furthermore,  if  we  let 
then  for  constant  U)  = ou  from  (2.22)  - (2.24),  we  have 


(2.23) 


(2.24) 


6 = A/vx  , 


(2.25) 


for  the  perturbation  in  wavenumber. 

We  have  discussed  first  order  perturbation  theory  only  because  that  is  all 
we  are  interested  in  here.  For  a discussion  of  second  and  higher  order  perturba^ 
tion  theory  see  Ref. 9. 

Although  the  perturbation  theory  developed  here  is  only  for  purely  real 
perturbations  of  real  eigenvalues,  it  may  readily  be  formally  generalized  to 
the  case  of  purely  imaginary  perturbations  of  real  eigenvalues.  This  formal 


generalization  is  accomplished  simply  by  requiring  the  T„,  u_.,  and  cp  to  be 


imaginary  while  the  T^.,  g^.',  D.  and  f^  are  real  in  Eq.  (2.18).  Then  the  H in 

ID’  3 1 M- 


(2.18)  is  purely  imaginary  and  the  i is  factored  out  before  the  integration  is 
performed. 


8. 


3 . Attenuation  Due  to  Air  Loading 

When  the  surface  of  the  substrate  is  in  contact  with  a fluid  medium  energy 

is  transmitted  from  the  acoustic  surface  wave  to  the  fluid.  Since  the  fluid  has 

a much  smaller  mechanical  impedance  than  the  solid,  the  resulting  attenuation  can 

readily  be  calculated  by  finding  the  steady  wave  induced  in  the  fluid  by  the 

small  displacement  at  the  surface  of  the  substrate  due  to  the  known  acoustic 

surface  wave  propagating  along  a free  surface  and  employing  (2.25)  with  (2.21) 

and  (2.18).  Since  the  fluid  experiences  small  wave  motion  about  a rest  position, 

we  may  employ  the  linear  equations  for  the  small  motion  of  a viscous,  compressible 
19 

fluid  , which  may  be  written  in  the  form 


where 


T.  . . = p u . , 

11.  i an’ 


a a c 2 a.  .a  a • a * a 

T . . = a u.  , 6 . . - — p 6 . . u + p (u.  ,+u.  .) 

n k,k  3 n k.k  ^ i.i  ’ 


il  V 


L>  3 1, 


and  A.3,  p3  and  p are  the  modulus  of  compression,  coefficient  of  viscosity  and 

cl 

mass  density  of  the  fluid,  u3  denotes  the  components  of  the  mechanical  displace- 
ment vector  of  the  fluid  and  6. . is  the  Kronecker  delta. 

il 

A solution  satisfying  the  differential  equations  and  boundary  conditions 
of  linear  piezoelectricity  for  acoustic  surface  waves^ ’ may  be  written 


in  the  form 


v-  , \ , n ^ i§  x -out 

\ „ (m)  (m)  m2  1 

u . = ) C A . e e 

3 L,  3 

m=l 

i,  , , , , iP  5x  i§  (x  -u>t) 

V (m)—  (m)  m2  1 

9 - ^ C A4  e e > 


for  propagation  in  the  direction  x^  with  x^  normal  to  the  surface  as  shown  in 
Fig.l.  For  a given  set  of  boundary  conditions,  say  traction-free  and  either 
short  circuit  or  open  circuit  electrical  conditions,  values  of  A_jm', 

A^m^  and  3 are  determined  numerically^ ’ 14’ 20  22.  These  calculations  have  been 


9. 


performed  for  various  cuts  of  lithium  niobate^' 22  and  quartz"0, 22  including 
22  24 

ST-cut  quartz  ’ and  are  assumed  known  for  purposes  of  this  work. 

Since  the  amplitude  of  the  mechanical  displacement  u3  in  the  substrate  is 
negligible  compared  to  u^  and  u^,  we  ignore  u^  in  the  treatment  of  the  fluid 

and  take  the  solution  in  the  fluid  in  the  form 

2 


a v „<»>,»  °^X2  i(?xru,t) 

U!  = l B ^ e 
m=l 

V-  , ' , S a 5x_  i(§x  -out) 
a \ (m)  (m)  m2  1 


■ I 


"2  L B °2  G 

m=l 


(3.4) 


which  satisfy  (3.1),  with  (3.2),  provided 

iugj-^A.  - j iuun  + | iuup,  (iuuy,  + p^V2)  + (x  - j iuup  - Pav2)(^  “ ^ ~) 

2 

- iuu  ^ ) cr"  + A.  - 4 iujp.  - p V2  [ ioup  + p v2]  = 0 , 

3/  J m L 3 3_J  3 


(3.5) 


and 


where 


(m)/  (m)  _ (X  - ioupt/3 ) icy  /[A  - (4icup./3)  + ioup,  a2  - p V2]  , (3.6) 

13  m m a 


v = u)/§  . 


(3.7) 


For  a given  uu  and  surface  wave  velocity  V known  from  the  aforementioned  surface 
6,14,20-24 


Eq.  (3.5)  is  a biquadratic  in  a and  yields  four  complex 

m 


wave  solutions 

roots.  The  two  corresponding  to  propagation  away  from  the  substrate  in  the 


x^-direction  are  inserted  in  (3.4)  and  (3.6)  for  the  solution.  Moreover,  these 
two  solution  functions  decay  with  respect  to  distance  from  the  surface  of  the 
substrate  in  the  x0-direction.  The  constants  are  determined  from  the  values 

of  u^  and  u.,  at  the  surface  of  the  substrate,  which  is  at  x^  = 0.  For  the 
solution  in  (3.4),  the  nontrivial  components  of  the  traction  vector  on  planes 


a . , a , a , 

T 21  ~ “ (U2, 1 Ul,  2* 


normal  to  x,  are  given  by 


10. 


= (A.  + 2i(u p,/3 ) u^  ^ + (A  - 4i(JU|i/3)  u^9^  • 


22 


(3.8) 


In  applying  Eq.  (2.18)  and  evaluating  the  normalization  factor  N in  (2.12) 

CM' ) 

the  integration  is  from  the  surface  of  the  substrate  at  x2  = 0 to  infinity  with 
respect  to  depth  and  over  a wavelength  of  the  acoustic  surface  wave.  The 
normalized  acoustic  surface  wave  eigensolution  is  given  by 


gj  = Uj/Nl  ’ ^ = cP/Ni  » 

where  u_.  and  cp  are  given  in  (3.3),  and  from  (2.12)  and  (3.3)  we  have 


(3.9) 


2 . tt  r 

Ni  = 1 3 p l 


4 4 C(m)A.(m)C(n)\(n)* 

k k 


L 

m=l  n=l 


(0  - p*  ) 

m n 


(3.10) 


and  represents  complex  conjugate.  The  first  perturbation  of  the  eigenvalue 
may  now  be  calculated  from  (2.18),  which  takes  the  form 

2TT/5 

f a 1 a I-! 

(3.11) 


Hl=- 


■ ' -J 

I [ 

o x2=0 


T2igi  + T22g2 


dx. 


because  T., ^ = 0 at  x2  = 0 in  the  eigensolution  in  which  the  correct  electrical 
continuity  conditions  have  been  employed  so  that  the  electrical  terms  in  (2.18) 
cancel  and  the  volumetric  perturbation  terms  T^  and  vanish  identically  in 
this  problem.  In  accordance  with  the  discussion  at  the  end  of  Section  2,  wr  now 
substitute  the  real  part  of  (3.9)  and  the  imaginary  part  of  (3.8)  into  (3.11), 

factor  out  an  i and  perform  the  integration  to  obtain 

2 


H = ilml  - tt  y B (-icon  D (m)Q'  + um  D y^g?  + (i  (A  + 2iou m- /3  ) n,( 
-L  L L^lm  z l 1 

(A  - 4i(ju|i/3)  D-!m)Q'm)92"|  ] , 


m=l 


x2=0 


(3.12) 
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where  the  differ  from  the  because  of  the  normalization  in  the 

normalized  surface  wave  solution  which  has  been  employed.  The  symbol  Im  denotes 
the  imaginary  part  of  the  complex  quantity  in  brackets  in  (3.12).  Note  that 
the  real  part  of  the  complex  quantity  in  brackets  in  (3.12)  gives  a small,  rela- 
tively unimportant,  change  in  surface  wave  velocity  which  we  do  not  bother  to 
present.  The  solution  we  have  obtained  can  be  improved  by  iterating  in  the 
manner  discussed  in  the  next  two  sections,  but  the  present  result  is  sufficiently 
accurate  that  it  is  not  purposeful  to  iterate  here. 

When  the  viscosity  p,  of  the  fluid  is  neglected  only  the  normal  components 
of  the  mechanical  displacement  and  traction  force  are  continuous  at  the  fric- 
tionless interface  and  in  this  case  the  solution  reduces  considerably.  First, 
there  is  no  sum  in  (3.4)  because  when  p = 0 (3.5)  yields  only  one  root  corre- 
sponding to  a wave  propagating  away  from  the  substrate.  Then  one  need  only 
eliminate  the  sum  in  (3.12)  and  set  = 0 wherever  it  occurs.  Under  these 
circumstances  reduces  to  the  considerably  simpler  form 

Hj  = irr\(aj2/§2)g2g*/[  (u)2/§2)  - (X/p3)]*5  . (3.13) 

Calculations  have  been  performed  for  6 as  a function  of  frequency  using 
both  (3.12)  and  (3.13)  and  the  results  are  plotted  in  Fig. 2.  As  can  be  seen  from 
the  figure,  with  increasing  frequency  there  is  increasing  deviation.  At  a fre- 
quency of  1GHz,  the  attenuation  calculated  from  (3.12)  is  in  very  close  agreement 

with  the  measurements  of  Szabo  and  Slobodnik'*'^  and,  of  course,  the  attenuation 

25 

calculated  from  (3.13)  is  not.  Bechmann's  constants  for  quartz  were  used  in 

the  calculation.  The  constants  employed  for  the  air  at  standard  temperature  and 

26 

pressure  (20 °C  and  1 atmosphere)  are 


p = 1.21  H X = 14.22  X 104  iL 
3 M M2 


18.22  X 


IQ’6  ™ 
M 


I. 


(3.14) 


12. 


Similar  calculations  have  been  performed  for  lithium  niobate,  and  the  results 

are  plotted  in  Fig. 3,  which  shows  reasonable  agreement  with  the  measurements  of 

5 

Slobodnik  . 


4 . Velocity  Dispersion  Due  to  Thin  Films 

When  a thin  film  of  a different  material  is  plated  on  a substrate  the  funda- 
mental essentially  Rayleigh  type  surface  wave  becomes  dispersive.  Moreover, 

14  27 

when  the  film  is  very  thin  compared  with  a wavelength,  it  has  been  shown  ’ 

28 

that  approximate  thin  plate  equations  can  be  used  w'th  great  accuracy  to 

treat  the  entire  mechanical  effect  of  the  plating  material  as  a boundary  condi- 

29  7 

tion  at  the  surface  of  the  substrate.  Both  Skeie  and  Auld  have  used  pertur- 
bation equations  to  calculate  the  initial  slope  of  the  Rayleigh  wave  dispersion 
curve  due  to  a thin  film.  In  this  section  we  essentially  repeat  this  calcula- 
tion in  order  to  emphasize  the  extreme  accuracy  with  which  the  initial  slope 
of  the  dispersion  curve  due  to  a thin  surface  film  can  be  determined  by  means 
of  perturbation  theory.  However,  we  go  one  step  further  and  show  how  the  equa- 
tion for  the  first  perturbation  of  the  eigenvalue  can  be  used  iteratively  to 
calculate  with  great  accuracy  considerably  more  than  the  initial  slope. 

The  aforementioned  approximate  plate  equations  for  a thin  isotropic  film 
plated  on  a substrate  may  be  written  in  the  form 


T . = - 6 . 2h  V 1 j--  )u  + u 1 + 2h,p/u.  , 

2d  Db  L\\'+2)jl  • a,ab  b,  aa  J r D ’ 


(4.1) 


where  p 2h'  and  H ' are  the  mass  density,  thickness  and  Lame  constants, 
respectively,  for  the  film  material,  the  subscripts  a and  b can  take  the 
values  1 and  3 but  not  2 and  denotes  the  components  of  the  traction  vector 
at  the  surface  of  the  substrate,  which  may  be  obtained  from  (2.3). 
application  Eq.  (2.18)  takes  the  form 


In  this 
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2TT/S 

Hl  = - I [T2j^j]dXl- 


(4.2) 


° V° 


e 6 

because  the  volumetric  terms  T\  ^ and  vanish  identically,  the  boundary  condi- 

1 - 1 

tions  employed  xn  the  surface  wave  solutions  are  = 0 and  f =0  and  cp  = 0 

in  the  conducting  thin  film.  Substituting  the  real  parts  of  (3.3)  for  the 
electrically  shorted  case  and  (4.1)  into  (4.2),  taking  the  u_.  in  (4.1)  to  be 
the  normalized  u_.  of  the  surface  wave  solution  in  (3.3)  evaluated  at  the  surface 
of  the  substrate  and  performing  the  integration,  we  obtain 

M ,a)2'\  » p'ui"  * / , p 

, , — + iM-  - - 

(*'  +2M.')  f2  > 1 1 5' 

where 

V = 2h'§  , 


V( . / (X ‘ + n ' ) p ,a)2'\  » p ;cu2  # [ , pV’-1) 

1 = " ~ ~Jgigl  - — 92g2  + VX  )g3g3J'  (4-31 


(4.4) 


and  the  g_.  are  the  components  of  the  normalized  displacement  vector  u^  at  the 
interface.  The  initial  slope  of  the  phase  velocity  dispersion  curve  for  an 
aluminum  film  on  ST-cut  quartz  has  been  calculated  from  (4.3)  and  is  plotted 
in  Fig. 4 along  with  the  correct  curve.  it  is  clear  from  the  figure  that  the 
agreement  of  the  initial  slopes  is  excellent. 

If  the  phase  velocity  determined  from  (4.3),  with  (2.21)-  (2.23).  for  a 
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small  nonzero  value  of  \ is  inserted  in  the  solution  for  ST-cut  quartz  ’ , 

with  the  traction  boundary  conditions  in  (4.1)  replacing  the  traction  free 
conditions  employed  in  the  treatments,  the  system  resulting  from  the  differ- 
ential equations  may  readily  be  solved  for  the  four  (3^  and  the  accompanying 
amplitude  ratios  and  then  the  system  resulting  from  any  two  of  the  three 

traction  boundary  conditions  along  with  the  condition  of  vanishing  potential 
may  readily  be  solved  for  the  amplitude  ratios  The 


resulting  solution  may  then  be  employed  in  (2.18),  which  now  takes  the  form 


14. 


H 


1 


2n/g 

j LT2j^  - T2 jU j]  «1  • 

° x =0 


(4.5) 


to  calculate  a value  of  phase  velocity  at  a nonzero  value  of  y.  Equation  (4.5) 
may  then  be  used  iteratively  to  calculate  the  phase  velocity  at  a given  value 
of  y to  any  desired  degree  of  accuracy,  subject  only  to  the  limitations  inherent 
in  (4.1).  To  this  end  we  substitute  the  real  parts  of  (3.3)  for  the  electric- 
ally shorted  case,  (4.1)  for  T . and 


1 E 1 .1 

T2j  = C2jUgk ,1  + ek2jf,k  ’ 


(4.6) 


into  (4.5),  take  u_.  as  before  and  perform  the  integration  to  obtain 

( / , U)2\  11*  tt  1*  * 1 , 

H1  Tf  Y -p  ^2/y35 * * * 93  + IF  (U3T23  +U3T23>’ 


(4.7) 


where  it  has  been  assumed  that  the  one  boundary  condition  that  has  not  been 
satisfied  exactly  is  the  condition  on  the  traction  T For  nonzero  y Eg.  (4.7) 
may  now  be  used  iteratively  to  calculate  the  phase  velocity  to  any  desired 
degree  of  accuracy.  The  phase  velocity  dispersion  curve  for  an  aluminum  film 
on  ST-cut  quartz  has  been  calculated  in  this  manner,  and  the  resulting  i urve  is 
indistinguishable  from  the  correct  curve  shown  in  Fig. 4.  Similarly,  the  phase 
velocity  dispersion  curve  for  an  aluminum  film  on  Y-Z  lithium  niobate  has  been 
calculated  and  the  results  are  plotted  in  Fig. 5. 


5 . Attenuation  Due  to  the  Electrical  Conductivity  of  Thin  Films 

When  a thin  conducting  film  is  plated  on  the  surface  of  a piezoelectric 

substrate  carrying  an  acoustic  surface  wave,  the  attendant  electric  field  causes 

currents  to  flow  in  the  film.  Since  the  film  is  thin  compared  with  a wavelength, 
the  full  electrical  conduction  equations  and  associated  boundary  conditions  need 

not  be  satisfied  for  the  film  and  an  approximate  thin  plate  conductivity 


15. 


equation  can  be  obtained  as  noted  in  the  Introduction.  If  we  assume  that  the 

electric  potential  essentially  does  not  vary  across  the  thickness  of  the  thin 

isotropic  film  and  that  the  electric  field  E.  is  derivable  from  the  electric 

1 

potential  cp  in  the  usual  way,  i.e.,  by 

Ei  = -V’  (5>1) 

then  application  of  the  conservation  of  charge  to  the  element  of  plate  shown 
in  Fig. 6 yields 

2h ' j>  NflJa  ds  + J D2  dS  = 0 , (5.2) 

C A 


where  N denotes  the  components  of  the  unit  normal  in  the  plane  and  J denotes 
a a 

the  components  of  the  current  in  the  plane.  Substituting  from  Ohm's  law 


Ja  = aEa  > 


(5.3) 


(5.1)  and  employing  the  divergence  theorem  in  the  plane,  we  obtain 


^h'acp^  +D2  = 0, 

where  a is  the  electrical  conductivity  and  D.,  may  be  determined  from 
and  the  surface  wave  solution  functions  in  (3.3).  Substituting  from 
into  (5.4),  we  obtain 


(5.4) 


(3.3) 


cp  (0 ) = iU)D2  (0)/2h'a£2  , (5.5) 

where  the  notation  (0)  means  evaluated  at  x2  = 0,  and  the  quantity  2h  'a  is 
called  the  sheet  conductivity  of  the  thin  surface  film. 

Even  though  we  can  readily  include  the  mechanical  influence  of  the  film 
by  means  of  Eq.  (4.1),  in  this  section  we  ignore  this  mechanical  effect  for 
reasons  of  expediency.  Under  these  circumstances  Eq.  (2.18)  takes  the  form 
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Hl=- 


o x2=0 


0 0 

because  the  volumetric  terms  T„  and  vanish  identically  and  the  mechanical 
terms  have  been  ignored.  At  this  point  it  should  be  noted  that  had  the  mech- 
anical terms  been  included,  Eq.  (5.6)  would  have  resulted  anyway  if  the  correct 
dispersive  surface  wave  solution  were  employed  because  then  (4.1)  would  be 
used  for  both  T . and  . in  (2.18).  However,  we  do  not  bother  with  this 
refinement  here.  For  future  use  we  note  that  the  complex  equation  (5.5)  yields 
the  purely  real  and  purely  imaginary  equations 

cp(0)  + cp*(0)  = — — - (D  (0)  - d!(0))  , (5.7) 

2h  'ar 

9(0)  - 9*  (0)  = r (0)  + D (0) ) , (5.8) 

2h,or  2 2 

which  we  will  have  occasion  to  use  separately  in  the  sequel. 

Since  we  are  interested  in  obtaining  results  over  an  extremely  large  range 
of  values  of  sheet  conductivity  2h/a,  we  must  perturb  from  the  short  circuit 
solution  for  large  values  of  sheet  conductivity  and  the  open  circuit  solution 
for  small  values  of  sheet  conductivity.  To  this  end  we  first  perturb  from  the 
short  circuit  surface  wave  eigensolution,  in  which  f^  = 0,  and  in  accordance 
with  the  discussion  at  the  end  of  Section  2,  consider  a purely  imaginary  H^. 

in  which  case  (5.6)  takes  the  form 

2TT/C 

i r * -xi  -ii* 

H1  = 4 .1  (9  - 9 ) 02  + ) dxi  • (5.9) 

O 


Substituting  from  (5.8)  into  (5.9)  and  performing  the  integration,  we  obtain 


Hj.  = (i0UTT/2h/CTf3)i2(O)32'f(O) 


(5.10) 
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where 


U _ 1 „S  U 

^2  ~ e2kXgk,£  ~ 2kf, k 


(5. 11) 


The  attenuation  calculated  from  (5.10)  at  a frequency  of  100  MHz  for  ST-cut 
quartz  is  plotted  as  the  curve  marked  short-circuit  in  Fig.  7.  This  curve  is 
valid  for  extremely  large  values  of  sheet  conductivity  only,  which  happens  to 
be  the  region  of  practical  interest.  For  somewhat  smaller  values  of  sheet  con- 
ductivity, a considerably  more  accurate  value  of  attenuation  may  be  obtained  in 
the  following  manner.  First  the  corrected  eigenvalue,  including  the  determined 
time  attenuation  as  well  as  the  phase  velocity,  may  be  substituted  in  the  system 
resulting  from  the  solution  of  the  piezoelectric  differential  equations,  which 
may  then  readily  be  solved  for  the  four  and  the  accompanying  amplitude 
ratios  A^m^.  Then  the  linear  algebraic  equations  resulting  from  the  three 
traction  boundary  conditions  may  be  solved  for  the  amplitude  ratios  C ' :C  ' : 

C<3)  : C^.  The  resulting  solution  may  then  be  employed  in  (5.6).  which  for  a 
purely  imaginary  H*  we  take  in  the  form 


i i r r * *i  'll*  * -'i*  ~i 

H1  = 4 J [(9-9  >tf2+£,  > ~ (D2-°2)(f  +f  >Jcbci  , 

o 


(5.12) 


in  order  that  the  iterative  perturbation  procedure  tend  to  satisfy  the  boundary 
condition  (5.5),  i.e.,  in  order  that  h[  vanish  for  the  exact  solution  that 

satisfies  (5.5).  For  a purely  real  H^,  (5.6)  takes  the  form 

2n/§ 

H1  = T J [(9+9*)  (^2+^')  " (D2  + D2)  (f1  + fly>]  dxi  • (5-13) 

o 

Substituting  from  (5.8)  into  (5.12)  and  (5.7)  into  (5.13)  and  performing  the 


integrations,  we  obtain 
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I _ iOKT  *111*  ^ ,tl*  n 

H1  T~ 77  ^2  ~ 2%  (f  °2  " °2f  * ’ 

2h  CT§  5 

R TT  A1  * ~ l* 

H1  = " 2?  (f  °2  + V > > 


(5. 14) 


(5.15) 


where  D.,  is  to  be  taken  to  be  the  same  as  & given  in  (5.11).  Equation  (5.14) 
determines  the  change  in  attenuation  from  any  given  state  and  Eq.  (5.15)  deter- 
mines the  change  in  phase  velocity  or  frequency  from  any  given  state.  Equa- 
tions (5.14)  and  (5.15)  may  now  be  used  iteratively  to  calculate  the  attenuation 


and  phase  velocity  to  any  desired  degree  of  accuracy  as  long  as  the  attenuation 
is  not  too  large.  The  attenuation  calculated  in  the  foregoing  manner  at  a 
frequency  of  100  MHz  for  ST-cut  quartz  yields  the  curve  below  the  short-circuit 
curve  in  Fig. 7.  The  iterated  attenuation  calculation  does  not  converge  for 
values  of  sheet  conductivity  lower  than  the  value  at  which  the  solid  curve  ends. 
We  believe  that  this  difficulty  arises  because  the  perturbation  procedure  wo 
employ  is  for  perturbations  off  real  eigenvalues  only,  and  in  the  region  in 
which  we  caraiot  perform  the  calculation  the  attenuation  is  too  large  for  the 
eigenvalue  to'  be  treated  as  approximately  real  for  the  accuracy  required. 

We  now  perturb  from  the  open-circuit  surface  wave  eigensolution,  in  which 


■1  C-1 

'2  ~ *Vf  ’ 


(5. 16) 


by  considering  a purely  imaginary  in  (5.6)  which,  with  (5.16)  takes  the  form 


2n/§ 

Hi  = i J - (d2 

o 


°2  > J L f +f  JdX! 


(5.17) 


Substituting  from  (5.7)  into  (5.17),  taking  cp  to  be  f and  performing  the  inte- 


gration, we  obtain 


H1  = (i§n2h'a/u))f1(0)f1'' (0)  . 


(5. 18) 


The  attenuation  calculated  from  (5.18)  at  a frequency  of  100  MHz  for  ST-cut 
quartz  is  plotted  as  the  curve  marked  open-circuit  in  Fig. 7.  This  curve  is 
valid  for  relatively  small  values  of  sheet  conductivity  only.  For  somewhat 
larger  values  of  sheet  conductivity,  a more  accurate  value  of  attenuation  is 
obtained  by  using  the  corrected  eigenvalue  in  the  governing  equations  in  the 
manner  set  forth  in  the  short-circuit  case  and  employing  the  resulting  solution 
in  (5.12)  and  (5.13),  respectively,  while  substituting  from  (5.7)  into  (5.12) 
and  (5.8)  into  (5.13)  and  performing  the  integrations  to  obtain 
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Equations  (5.19)  and  (5.20)  may  now  be  used  iteratively  to  calculate  the  atten- 
uation and  phase  velocity,  respectively,  to  any  desired  degree  of  accuracy  as 
long  as  the  attenuation  is  not  too  large.  The  attenuation  calculated  in  the 
foregoing  manner  at  a frequency  of  100  MHz  for  ST-cut  quartz  yields  the  curve 
below  the  open-circuit  curve  in  Fig. 7.  This  time  the  iterated  attenuation 
calculation  does  not  converge  for  values  of  sheet  conductivity  higher  than  the 
value  at  which  the  solid  curve  ends,  and  the  reasons  for  this  are  the  same  as 
those  discussed  earlier  in  the  short-circuit  case.  Similar  calculations  have 
been  performed  for  Y-Z  lithium  niobate  and  the  results  are  plotted  in  Fig. 8. 
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FIGURE  CAPTIONS 


Diagram  Showing  the  Free  Surface  of  a Semi-infinite  Solid 

Attenuation  of  Surface  Waves  on  ST-Cut  Quartz  due  to  Air 
Loading  at  Standard  Temperature  and  Pressure  (a)  with  and 
(b)  without  Fluid  Viscosity.  The  dot  is  the  value  measured 
by  Szabo  and  Slobodnik 

Attenuation  of  Surface  Waves  on  Y-Z  Lithium  Niobate  due  to 
Air  Loading  at  Standard  Temperature  and  Pressure  (a)  with  and 
(b)  without  Fluid  Viscosity.  The  dot  is  the  value  measured 
by  Slobodnik5 . 

Lowest  Straight-Crested  Phase  Velocity  Dispersion  Curve  for 
an  Aluminum  Film  on  ST-Cut  Quartz 

Lowest  Straight-Crested  Phase  Velocity  Dispersion  Curve  for 
an  Aluminum  Film  on  Y-Z  Lithium  Niobate 

Diagram  Showing  an  Arbitrary  Element  of  the  Thin  Conducting 
Film 

Attenuation  of  Surface  Waves  on  ST-Cut  Quartz  due  to  Thin 
Conducting  Films  as  a Function  of  Sheet  Conductivity  at  a 
Frequency  of  100  MHz 

Attenuation  of  Surface  Waves  on  Y-Z  Lithium  Niobate  due  to 
Thin  Conducting  Films  as  a Function  of  Sheet  Conductivity  at 
a Frequency  of  100  MHz 
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